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ABSTRACT 

Wc propose a mechanism for the superhump modulations observed in optical pho- 
tometry of at least two black hole X-ray transients (SXTs). As in extreme mass-ratio 
cataclysmic variables (CVs), superhumps are assumed to result from the presence 
of the 3:1 orbital resonance in the accretion disc. This causes the disc to become 
non-axisymmetric and precess. However the mechanism for superhump luminosity 
variations in low mass X-ray binaries (LMXBs) must differ from that in CVs, where it 
is attributed to a tidally-driven modulation of the disc's viscous dissipation, varying 
on the beat between the orbital and disc precession period. By contrast in LMXBs, 
tidal dissipation in the outer accretion disc is negligible: the optical emission is over- 
whelming dominated by reprocessing of intercepted central X-rays. Thus a different 
origin for the superhump modulation is required. Recent observations and numerical 
simulations indicate that in an extreme mass-ratio system the disc area changes on 
the superhump period. We deduce that the superhumps observed in SXTs arise from 
a modulation of the reprocessed flux by the changing area. Therefore, unlike the situa- 
tion in CVs, where the superhump amplitude is inclination-independent, superhumps 
should be best seen in low-inclination LMXBs, whereas an orbital modulation from 
the heated face of the secondary star should be more prominent at high inclinations. 
Modulation at the disc precession period (10s of days) may indicate disc asymmetries 
such as warping. We comment on the orbital period determinations of LMXBs, and 
the possibility and significance of possible permanent superhump LMXBs. 

Key words: accretion, accretion discs - stars: binaries close - stars: individual: 4U 
1915-05, GX 9+9, GS 1124-68, GRO J0422-f32, GS 2000+25, GRS 1716-249 



1 INTRODUCTION 

Superhumps are periodic optical modulations observed in 
superoutbursts of the SU UMa dwarf novae (see Warner, 
1995 for a review). Their most striking property is that 
their periods Psh are slightly longer than the orbital period 
PoTh, typically by 1 — 7 %. Since the work of Whitehurst 
(1988), Whitehurst & King (1991) and Lubow (1991a, b) 
it is now understood that superhumps are a consequence 
of the presence of the 3:1 orbital resonance within the ac- 
cretion disc. This causes the disc to become eccentric, and 
to undergo slow prograde precession in the inertial frame. 
The secondary star thus repeats its motion with respect to 
the disc on the beat period between the orbit and this pre- 
cession, which is therefore slightly longer than Porb- This 
relative motion of the secondary star modulates the disc's 
streamline geometry on the superhump period, causing its 
viscous dissipation to vary on the same period. This intrinsic 
variation of the disc light satisfies one of the basic observa- 



tional features of CV superhumps, namely that their oc- 
currence is independent of binary inclination (Warner 1995, 
section 3.6.4.1). 

The defining resonance condition restricts the mass ra- 
tio q — M2/M1 to extreme values 



q <, 0.33. 



(1) 



Since in most CVs the white dwarf mass Mi lies in a narrow 
range (~ 0.6 — O.8M0), and M2 is often strongly correlated 
with the orbital period P = Phr hr (i.e. M2/M0 ~ O.llPhr), 
this confines the occurence of superhumps to short orbital 
periods, mostly below the CV period gap (so P^r ^ 2). In 
practice almost all of these systems are dwarf novae, with 
the superhumps occurring during superoutbursts. However 
persistent systems satisfying the resonance condition (|l|) do 
exist, and show the superhump phenomenon permanently 
('permanent superhumpers', Patterson 1999). 'Negative su- 
perhumps', in which Psh < Porb, also exist, and clearly cor- 
respond to retrograde disc precession in the inertial frame. 
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There is no currently accepted explanation for the retro- 
grade precession inferred from negative superhumps in CVs. 

Murray (2000) shows that the precession period, and 
thus the quantity 



(2) 



depend on both the mass ratio and other conditions such 
as the disc pressure or temperature. Small mass ratios q 
produce small e, as does a high disc temperature. Taking 
this to extremes, an extreme mass ratio system with a very 
hot disc could have a disc with retrograde precession. 



2 SUPERHUMPS IN LMXBS 

The resonance condition (1) is easily fulfilled in binaries with 
large primary masses Mi . The prevalence of superhumps in 
outbursting CVs suggests that soft X-ray transients, par- 
ticularly those containing black holes, (see Marsh 1998 for 
a review) might be prime candidates for showing super- 
humps. Recognising this, O'Donoghue and Charles (1996) 
carefully surveyed the available observational evidence, and 
concluded that superhumps had been seen in outbursts of at 
least two black-hole SXTs (GS1124-68 = X-Ray Nova Mus 
1991, GRO J0422-h32, and probably GS 20004-25). At about 
the same time, superhumps were claimed in observations 
from the 1993 outburst of GRS 1716-249 (=X-Ray Nova Oph 
1993 = GRO J1719-24, Masetti et al 1996), but as the orbital 
period of this system is unknown, this claim cannot be rig- 
orously assessed. It is notable that the superhump excesses 
e are only ~ 1 — 2%, which with Porb ?5 5hr, implies much 
longer disc precession periods, Pprcc = Psh/e ~ 10 — 50 days, 
than in CVs. 

Although supporting the importance of the 3:1 orbital 
resonance, the existence of superhumps in LMXBs presents 
us with an apparent difflculty. The intrinsic dissipation in 
LMXB discs has long been known to be a negligible contrib- 
utor to their optical light (e.g. van Paradijs and McClintock 
1995). The argument is simple: the X-ray luminosity Lx 
gives an estimate of the accretion luminosity r)Mc^ (where 
rj ~ lO'^" erg g~^ is the efficiency of rest-mass conversion) 
and thus the accretion rate M on to the central star (neu- 
tron star or black hole) in an LMXB. This immediately 
gives an estimate of the optical luminosity I/opt(visc) of a 
steady-state accretion disc surrounding the central object 
(cf Frank et al., 1992). For all persistent LMXBs the ratio 
iopt(visc)/Lx predicted by this method is far smaller than 
the observed value (van Paradijs & McClintock 1995). This 
conclusion can be extended to SXTs in outburst, as their 
discs can then be regarded as approximately steady. 

The explanation for the excess optical luminosity of 
LMXBs is straightforward. If the disc intercepts even a small 
fraction of the central X-ray luminosity, this completely 
dominates its own intrinsic dissipation. For a point source 
at the centre of the disc, the irradiation temperature Tin- is 
given by 



4 _ vMc'{l-f3) H 



(e.g. van Paradijs, 1996) Here /3 is the X-ray albedo, H{R) 
is the local disc scaleheight, and 



dlnH , 
^=1 dhTR-' 



(4) 



Viscous dissipation alone gives an effective temperature 
Tvisc, with (e.g. Frank et al. 1992) 



3GMM 



at disc radii R much larger than the radius, r« , of the central 
object (mass Al). Dividing, one finds 



rp4: 



H 



(6) 



For LMXBs, the combination (27)c^/3GM)(l - fi) is of or- 
der rZ^ , where r, is the radius (e.g. event horizon) of the 
central star. Thus for a large enough disc, i.e. one with 
R » rt{R/ H g), irradiation wins over intrinsic dissipation 
despite the small solid angle ~ H/R of the disc, because 
Tj^j falls off only as R~^ , whereas T^i^^c goes as P"^ This 
condition is very easily satisfied in all LMXBs {R/r, 10'*, 
while R/Hg ^ 10^). Lopt is thus predominantly a result of 
disc irradiation. In agreement with this, van Paradijs & Mc- 
Clintock (1994) show that for a sample of 18 LMXBs the 
observed Lopt scales with Lx and disc size approximately 
as expected. 

It is perhaps worth noting that the mere existence of 
efflciently irradiated discs is a challenge to theory, which 
usually predicts that a disc heated by central X-rays adopts 
a convex shape (formally g < 0), and thus shields most of its 
area from the central flux (e.g. Cannizzo 1994, Dubus et al 
1999). However the observational evidence that LMXB discs 
are irradiated is overwhelming, and we adopt this view here 
(cf King & Ritter, 1998). 

The overwhelming dominance of irradiation over in- 
trinsic dissipation means that the explanation of the su- 
perhump luminosity variations in CVs will not work for 
LMXBs. The intrinsic dissipation at a given disc radius R, 
e.g. the resonant radius, varies as the local accretion lu- 
minosity I/acc(P) = GMM/R: the compactness M/r^ of 
the central object is irrelevant. The intrinsic superhump lu- 
minosities Lsh in LMXBs and CVs are thus in the ratio 
(M/P)lmxb/(M/P)cv. Even for black-hole systems, where 
Mlmxb/Mcv ~ 10, the longer orbital periods and thus 
larger values of R make this ratio of order unity. We conclude 
that intrinsic LMXB superhumps have luminosities Lsh sim- 
ilar to those in CVs. But the latter are a fraction / 0.1 
only of the total intrinsic optical disc luminosity Lopt (vise), 
and the same will be true in LMXBs. Thus intrinsic super- 
hump luminosity variations in LMXBs have amplitudes 



Lsh 
Lirr 



Lsh 



Lopt (vise) Li 



W(visc) ^ ^ ^ ^Q_3 ^ ^Q_4 



(7) 



of the observed disc brightness resulting from irradiation. 
Superhump variations powered by intrinsic viscous dissipa- 
tion are therefore negligible in LMXBs, and cannot explain 
the observed superhump amplitudes. 



3 DISC AREA VARIATIONS 

The work of the last Section shows that one can only res- 
cue the resonance theory of superhumps for LMXBs if the 
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Figure 1. Disc surface area as it varies over the course of a 
single superhump cycle. The solid line shows the change in surface 
area of the entire disc (+ stream), the dot-dash line shows disc 
region with densities > 1% of the maximum density, dotted line 
shows regions with density > 5% of the maximum and the dashed 
line shows regions with density > 10% of the maximum. In this 
simulation the superhump period is 1.0295Pojij . 



efficiency of tfie precessing disc in reprocessing the central 
X-rays varies on the superhump period. Writing 



(8) 



we see that this requires that the effective area A^s which 
the disc presents to the X-rays must vary on this cycle. In 
principle this could occur because the disc aspect ratio might 
vary on this cycle. Smale et al. (1992) suggest this occurs 
because the vertical component of gravity for the disc rim 
will be significantly increased when the outermost part of the 
elliptical disc coincides in azimuth with the mass donor star. 
They do not, however, present any quantitative assessment 
of how large an effect this will be. 

In contrast, by far the simplest possibility is that the 
total disc surface area varies. This may indeed have been 
seen in observations of a persistent superhumper, V348 Pup 
(Rolfe et al, 2000). Moreover, SPH simulations also predict 
a ~ 10% variation of the area. Figure ^ shows the disc 
area changes found in a simulation (Murray 2000) of the 
dwarf nova OY Car. This simulation assumed a mass ratio 
q = 0.102, similar to those measured in black hole SXTs, and 
found a superhump period excess e = 0.0295. The maxima 
of the area variations are in phase with the intrinsic vis- 
cous dissipation maximum, while the area minima lag the 
dissipation minimum by about 45 degrees. From (^ we see 
that the area variations of Fig. |l| give approximately the pre- 
dicted optical light curve, apart from the effect of dilution 
by other sources of optical emission. The most important of 
these is the X-ray heated face of the secondary star, whose 
effective temperature Tirr, 2 we can crudely estimate from 



irr, 2 



(9) 



with a, i?2 the binary separation and secondary radius. Since 
R2/a ~ 0.462[g/(l + q)Y^^ (Roche geometry), and the disc 
radius i? is a substantial fraction of a, we see from (^ that 
the heated secondary and the disc have comparable opti- 
cal brightness. In an LMXB observed at low inclination this 
effect will reduce the predicted fractional superhump ampli- 
tude somewhat below the variations seen in Fig. |l| At higher 
inclinations the changing aspect of the heated face will pro- 
duce an optical modulation at Porb which will compete with 
the superhump variation at P^h ■ Moreover the foreshortening 
of the disc at such inclinations will weaken the superhump 
modulation. We would therefore expect to see both modu- 
lations only in rare cases, and SXTs to divide into groups 
showing one or other of them. The similarity of and (^ 
suggests that these two groups should contain roughly equal 
numbers of systems, approximately as observed (see below). 
The neutron star transient XTE J2123-058 is at high or- 
bital inclination, and showed a clear optical modulation at 
Porb, comprised of the changing aspect of the heated face 
and grazing eclipses (Zurita et al 2000) . There was no trace 
of a superhump modulation, but this is to be expected for a 
neutron star system with Porb = 6.0 hr, since the mass ratio 
is unlikely to satisfy (1). 



4 DISCUSSION 

We have shown that variations of the disc surface area on 
the superhump cycle offer a plausible explanation for the su- 
perhump light curves observed in LMXBs. The superhump 
modulation in area (Figure 1) shows that the geometry of 
the accretion flow varies on the superhump period. There- 
fore: (i) the area of the disc visible to the observer changes, 
causing a modulation in the optical flux; (ii) the solid angle 
the disc subtends at the X-ray source will probably change, 
causing a modulation in the intercepted fraction of Lx . Both 
these factors contribute to the superhump modulation. 

The small observed values of e are a natural consequence 
of the extreme mass ratios q and high disc temperatures in 
LMXBs. The incidence of detected superhump variations, 
currently unambiguously seen in 2 or 3 SXTs, is in line with 
expectations from this model, as data able to address the 
issue was not collected for most SXTs. If the disc had any 
front-back asymmetry, perhaps as a result of warping under 
radiation-induced torques (Pringle 1996, Wijers and Pringle 
1999), detection of optical modulation on the disc precession 
period Pproc might be possible. However the small observed 
e values mean that Pproc is of order 10s of days, and thus fre- 
quently comparable with the duration of the outburst itself. 
Since neither Porb nor Pprcc will in general be known until 
the outburst has flnished, this highlights the importance of 
obtaining accurate photometry of any 'orbital' variations at 
as many stages of an outburst as possible. 

In SXTs we can establish Porb accurately from radial ve- 
locity measurements in quiescence, and thus confldently say 
whether a given photometric modulation is or is not a su- 
perhump. In a persistent LMXB, it is generally only possible 
to determine Porb if the system is at high orbital inclination 
so that X-ray eclipses or dips are seen as the central source 
is periodically occulted by the mass donor star or vertical 
structures in the disc rim. Of the five possible LMXB super- 
humpers listed in Ritter and Kolb (1998) four are SXTs. The 
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fifth is the dipping source V1405 Aql which was originally 
suggested as a possible superhumper by White (1989). We 
discuss this source and GX 9+9 individually, before giving a 
general discussion. 

4.1 V1405 Aql (A 1916-05, 4U 1915-05) 

V1405 Aql was comprehensively studied by Callanan, Grind- 
lay and Cool (1995, hereafter CGC). Their observational 
findings were 

(i) the optical period is 50.4589 mins 

(ii) the X-ray period is 50.00 mins 

(iii) in an "anomalously low optical state" (0.5 mag below 
normal) both 50.46 and 50.00 min periods were observed in 
the optical 

(iv) the 50.4589 min period seems to have phase stability 
over 7 years. 

CGC's conclusions were first that 50.4589 mins is the orbital 
period Porb, since the superhumps in outbursting SU UMa 
stars show neither period nor phase stability over weeks, let 
alone years. Second, in the anomalously low optical state 
the intrinsic (i.e. viscously generated) brightness of the disc 
should make a significant contribution, rendering "SU UMa 
- like precessing behaviour" more observable. Finally phase- 
wandering of the optical dips when folded on the X-ray pe- 
riod prevents simple interpretation of the periodicities. 

Since 1915-05 is a steady X-ray source, it is misleading 
to compare it with the SU UMa CVs in which the disc is 
subject to the thermal-tidal instability. SU UMa stars ex- 
hibit superhumps during their outbursts, when their discs 
are in the high viscosity state for days to weeks. SU UMa 
discs therefore change their radii and mass distributions on 
timescales of days, and their superhumps evolve as a con- 
sequence of this. For a steady disc satisfying equation (1) 
we should expect the superhumps to settle into a stable os- 
cillation: simulations of precessing eccentric discs settle into 
periodic behaviour unless the boundary conditions change 
(Murray 1998, Whitehurst private communication). We con- 
clude that the phase stability of the 50.4589 min period is 
no barrier to interpreting it as Psh_- We accordingly identify 
the 50.00 min period as Porb. 

The ratio of viscously-generated to irradiation- 
generated flux remains constant (equation |^ unless either 
the disc is gaining or losing mass at some radius, or the ge- 
ometry changes. In fact all the non-negligible components 
of the optical light curve: Lace, Lirr, and the irradiated 
flux from the mass donor, remain in the same ratios for 
a steady state accretion flow with any M unless the geome- 
try changes. Therefore it is likely that CGC's "anomalously 
low optical state" coincides with a change in the geometry 
of the accretion flow, making irradiation of the mass donor 
star more prominent. If the geometry did not change there 
would be no reason for the modulation at Porb to become 
more prominent (the ellipsoidal variation would introduce a 
signal at Porb/2, but is never likely to be observable as the 
intrinsic flux from the mass donor is certain to be negligible 
in a short period X-ray-emitting neutron star binary like 
V1405 Aql). 

Since we are proposing that the geometry of the disc 
changes on Pgh, the deep dips seen in figure 7 of CGC could 
well be analogous to the "superdip" seen in OY Car by 



Billington et al (1996). This interpretation can also explain 
the phase- wandering of the optical dips when they are folded 
on Porb = 50.00 mins, hence negating the third of CGC's 
conclusions. 

We identify V1405 Aql as a persistent irradiated super- 
humper. With our interpretation e = 0.009, and the inferred 
period of apsidal disc precession is Pprec = 3.8 days. Inter- 
estingly, the shape of the optical light curves was found to be 
modulated on a ~ 4 day period (Smale et al 1989). Homer 
et al (2000) have extensive X-ray and optical data which 
seem consistent with our interpretation and with warping of 
the disc. Nodal precession of a warped disc may introduce 
futher observable periodic modulations, analogous to the 35 
day period in Her X-1 (Scott & Leahy 1999), and may be 
responsible for the 199 day cycle in 4U 1915-05 reported by 
Smale (1994). 

4.2 GX 9-K9 (4U 1728-16, Oph X-1, 2S 1728-169) 

This system was suggested as a possible persistent su- 
perhumper by Haswell and Abbott (1994) who found a 
4. 1744 ±0.0002 hr modulation in I band photometry. Schae- 
fer (1990) had previously reported a B band modulation of 
period 4.198 ± 0.028 hr, and Hertz and Wood (1988) found 
an X-ray period of 4.19 ± 0.02 hr. These periods are clearly 
all mutually consistent. With Porb = 4.2 hrs and a neutron 
star primary, the system has a mass ratio q ^ 0.28 (Schaefer 
1990) satisfying (1), since the secondary mass M2 cannot 
exceed the mass ~ O.5M0 of a main-sequence star filling 
the Roche lobe (see below). We should therefore expect GX 
9+9 to harbour a persistently precessing elliptical disc. The 
~ 10% fractional amplitude of the optical modulations in 
GX 9+9 is similar to those reported by O'Donoghue and 
Charles (1996) for the superhumps in SXTs, and the changes 
in the shapes of the optical light curves shown in Haswell 
and Abbott (1994), figure 3, are suggestive of a light curve 
shape modulation at Pproc. 

Since there are extensive RXTE observations of this 
source in hand, a more precise X-ray period is likely to soon 
be available. Clearly this, and more extensive optical pho- 
tometry is required to test our hypothesis that GX 9+9 is 
an irradiated disc persistent superhumper. Kong et al (2000) 
will address this. 



4.3 Many LMXB superhumpers? 

Table 1 lists LMXBs in order of increasing Porb as reported 
in the literature. The third column in this table gives the 
nature of the modulation leading to the Porb determination. 
For the eclipsing systems and those in which orbital motion 
has been measured (denoted either "opt RV" for mass donor 
radial velocity modulations or "pulsation RV" for pulse tim- 
ing modulations) Porb is securely determined. The remaining 
17 period determinations arise from modulations in X-ray, 
UV, or optical fiux and cannot be identified as Porb with 
certainty; these 17 periods are in boldface. 

We can show that the condition (1) for superhumps is 
likely to hold for the first 11 systems listed in Table 1. The 
accretor presumably cannot be less massive than a neutron 
star, and indeed is known to be such a star in the first 10 
cases because of the presence of X-ray bursts. Accordingly 
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we can assume that Mi <; 1.4M0. The hkely candidates 
for the donors are either main-sequence stars, satisfying 
M2 — O.llPhrM© (e.g. King, 1988), or degenerate stars, 
obeying M2 ~ 0.015(1 + X^^'^P^^^^Mq (King, 1988) where 
X is the fractional hydrogen content by mass. (Degener- 
ate companions are hkely in the 5 systems in Table 1 with 
Phr < 1.) We thus find 

qMS ^ 0.079Phr, gdeg ^ 0.0115(1 + Xf/^p-\ (10) 

for the two cases, showing that indeed the condition (1) for 
superhumps holds for Porb ^ 4.2 hr. The only possible ex- 
ception to this would be a donor star which has undergone 
thermal-timescale mass transfer, as is thought to have oc- 
curred in for example Cygnus X-2 (King & Ritter, 1999). 
These stars are stripped down to their denser central re- 
gions, and have larger masses when filling the Roche lobe at 
a given period. However it is unlikely that periods 2 hr 
are accessible to this kind of evolution (King et al., 2000). 

Accordingly we expect that for the lower inclination 
systems with Porb ^ 4.2 hr, an optical or UV modulation 
is more likely to be a superhump than an orbital modula- 
tion. For high inclination systems satisfying (1) it is possible 
that the X-ray dip behaviour depends on the disc's preces- 
sional and superhump cycles, rather than straightforwardly 
indicating Porb. 

For a very extreme mass ratio, q < 0.02, the circular- 
ization radius, Pcirc, is bigger than the 3:1 resonance radius 
(c.f Frank et al. 1992, eqns (4.17, 5.75); Warner 1995 eqn 
(2.4a)). In this case we expect the disc to show persistent 
superhumps for all values of M. An LMXB with a lOM© 
black hole primary and a secondary ^ 0.2Mq would satisfy 
this, and have Porb — 2 hr. 

If it were possible to show unambiguously that an 
LMXB with Porb > 5 — 7 hr exhibited superhumps, this 
would suggest it harboured a black hole primary. We note, 
of course, that superhumps alone cannot provide proof of a 
black hole primary, as the secondary could have lower mass 
than a main sequence star, for example if it were slightly 
evolved. This might allow superhumps even though the pri- 
mary was a neutron star. However, coupled with a lack of 
Type 1 X-ray bursts, and with the presence of X-ray spectral 
and timing indicators of black hole candidacy, this line of ar- 
gument could prove useful, since it it extremely difficult to 
obtain a dynamical mass determination for the accretor in a 
persistent LMXB. It is notable that all the identified black 
holes among LMXB primaries (denoted BH in Table 1) are 
in transient systems. The nature of the compact objects in 
steady LMXBs which neither pulse nor burst is currently an 
undesirably open question. 
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